Adiabatic pumping in the quasi-one-dimensional triangle lattice 
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We analyze the properties of the quasi-one-dimensional triangle lattice emphasizing the occurrence 
of flat bands and band touching via the tuning of the lattice hopping parameters and on-site energies. 
The spectral properties of the infinite system will be compared with the transmission through a 
finite piece of the lattice with attached semi-infinite leads. Furthermore, we investigate the adiabatic 
pumping properties of such a system: depending on the transmission through the lattice, this results 
in nonzero integer charge transfers or transfers that increase linearly with the lattice size. 



I. INTRODUCTION 

Quantum pumping is a process, where a periodic vari- 
ation of system parameters leads to a density flux of par- 
ticles through the system in spite of the absence of an ex- 
ternal voltage. The first proposal for quantum pumping 
was formulated by Thouless in 1983, 1 and it considered 
particles in a periodic system described by Bloch wave- 
functions. The Thouless approach is restricted to insulat- 
ing systems and here the number of particles transferred 
in one period is always quantized to integer values. 

Another pumping approach is based on elastic scat- 
tering matrices. It was formulated by Brouwer 2 and 
Biittiker. 4 6 and allowed to explain experimental results 
by Switkes 3 for finite open systems. In general, the par- 
ticle transfer within a pumping cycle is found not to be 
quantized to integer values. Both approaches are equiv- 
alent if transmission through the finite static system is 
suppressed along the whole pumping cycle 13 and in this 
case, the charge transfer is quantized. 

In this Article, we investigate quantum pumping in 
a quasi-one-dimensional lattice, i.e. a one-dimensional 
lattice model with a basis of two (or several, in gen- 
eral) lattice sites. These kind of systems are basic 
model systems for the verification of fundamental phe- 
nomena such as the Aharonov-Bohm effect in presence of 
electron-interaction 7,8 , the Aharonov-Casher effect 9 , or 
the combination of the two effects. 10 In some recent ex- 
periment, the Aharonov-Casher effect has been observed 
in quasi-one-dimensional chains realized with topolog- 
ical insulators. 11 Quasi-one-dimensional linear chains 
have also been proposed for studying effects of spin- 
polarization in presence of spin-orbit interaction. 12 

Here, we will focus on the quasi-one-dimensional tri- 
angle lattice (Fig. 1). This lattice is characterized by 
a basis of two lattice sites. Its band structure exhibits 
flat bands and band touching as a function of the lattice 
parameters. First we work out the conditions for the ap- 
pearance of these peculiar features. Then, we relate the 
properties of the infinite lattice to the transmission prob- 
ability of a finite piece of lattice which is connected via 
semi-infinite leads to some particle reservoirs. Finally, 
we apply Brouwer's formalism 2 for adiabatic quantum 
pumping (AQP) to the scattering matrix of this finite-size 



structure. We find that the charge transfer is finite for 
a pumping parameter cycle surrounding a band touching 
configuration, even if along the cycle the transmission is 
inhibited. Further, if the pumping cycle traverses config- 
urations with nonzero transmission, the charge transfer 
has in average a linear dependency on TV, the number of 
unit-cells of the triangle lattice (see Fig. 1) — it increases 
linearly with the lattice length. However, this results is 
valid only for the charge transfer but not for the pumping 
current. In fact, in order to fulfill the adiabatic approxi- 
mation, the pumping frequency has to decrease with the 
number of unit-cells. 
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FIG. 1: Tight-binding triangle lattice with semi-infinite leads. 
A number of unit cells (N) with the two sites A and B are 
connected with hopping parameters t±, £2 and £3 to form the 
triangle lattice of length L — N a. 

II. THE TRIANGLE LATTICE 

A. Spectral properties 

We define the quasi-one-dimensional triangle lattice 
within the tight-binding formalism. It has a basis con- 
taining two lattice sites A and B with on-site energies ca 
and cb and hopping parameters £1, £2 and £3 as sketched 
in Fig. 1. For the infinite lattice — no leads — the en- 
ergy spectrum of the system can be obtained via Bloch's 
theorem. The Hamiltonian in reciprocal space reads 



ca + 2£ 3 cos k ti + t2e 



h + he 1 ' 



(1) 



with n = k a. The eigenvalues can be expressed with 
dimensionless variables as 

6 i,2(ft) = cos k ± ^ [Ae + T1T2 + cos k] 2 + r 2 + r| 
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FIG. 2: (Color online) Illustration of the band gap for the 
infinite triangle lattice. Brightness visualizes the value of the 
energy gap. White dashed curves correspond to the appear- 
ance of a flat band. The crosses (green) are located at band- 
touching (a) or band-crossing points (b). Corresponding pa- 
rameters in (a): T2 = 8; (b): Ae = 0.5. 



e i,2(ft) = cos k ± y [Ae + cos k] + r\ + r| + 2tiT2 cos ft 

" (2b) 

where all parameters are expressed in units of £3: = 
t</t 3 , e(«) = (^^)-(e A + e B )/2)/t 3 , Ae = (e A -e B )/2t 3 , 
we introduce here also the dimensionless coupling con- 
stant between the leads and the central system c = 
t^/(ti,ts). We consider the system in absence of external 
magnetic fields, hence, due to the conservation of time- 
reversal- symmetry all lattice parameters are real valued 
[H*(k) = H(k)]. However, we assume the possibility 
of sign changes of the hopping parameters and on-site- 
energies. Depending on the parameters, the energy gap 
between the two bands, which is denned as the differ- 
ence between the minimum of the upper and the maxi- 
mum of the lower band, is visualized via the color scale 
in Fig. 2(a)-(b). Here, the spectrum (2) exhibits a one- 
dimensional Dirac point [Fig. 3(a)], band crossing with 
fiat bands, and opening of gaps as visualized in Fig. 3(e) 
and 3(g), respectively. In order to observe a band touch- 
ing or crossing the square root in (2) must vanish for a 
specific n. This is realized for either of the three config- 
urations 



r\ = r 2 
r\ = -r 2 
Ti = r 2 = 



Ae = 1 
Ae = -1 
lAel < 1 . 



(3a) 
(3b) 
(3c) 



The first two cases are included in Fig. 2(a). On the other 
hand, case (3c) corresponds to the splitting of the lattice 
into a linear chain (A-sites) with a cosine band and the 



isolated B-sites leading to a flat band (see Fig. 3(e) and 
Fig. 2(b)). The parameter Ae is then responsible for the 
relative position between the flat band and the cosine 
band. If |Ae| < 1, the bands cross at two points within 
the first Brillouin zone. 

The general condition for the occurrence of a flat band 
is immediately obtained by (2a) and reads 



rjVf - 2rir 2 Ae = 



(4) 



and is visualized in Fig. 2(a)-(b) by the white dashed 
lines. 



B. Transport properties 

We consider a finite piece of the triangle lattice coupled 
to particle reservoirs via two semi-infinite one-channel 
leads, which are modeled via linear chains. The calcu- 
lation of the elastic scattering matrix can be done via 
the Fisher-Lee-relation 16 18 ' 24 — this relates transmis- 
sion and reflection amplitudes to the retarded equilib- 
rium Green's functions of the system. The elements of 
the Green's function of the infinite system required in or- 
der to determine the transport properties are determined 
through the calculation of the self-energies of the leads 
that are added to the lattice Hamiltonian. 17 The prob- 
lem is then reduced to a finite set of algebraic equations. 
The details of the calculation of the scattering matrix 
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FIG. 3: Collection of energy bands s(k) of the infinite lattice, 
transmission probability Tjv(e) and a(e) for the finite open 
lattice (N = number of unit-cells), (a)-(d): band touching 
with n = T2 = Ae — 1; (e),(f): Band crossing with n = T2 = 
0, Ae = 0.5 ; (g),(h): Splitting with n = r 2 = 0.1, Ae = 0.5. 
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are presented in App. A. The analytic expression for the 
5-matrix reads 



S12 — S21 

S 



-2i<7 



Bi 1 + (/?i+/3 2 )B 1 2-/?i/? 2 B 2 2 
2icr(B 12 -ffjBm) 



(5a) 



Bn + (/3i + /3 2 )B 12 -/3 1 /3 2 B 22 



1 . (5b) 



with j G {1,2} are the left and the right lead, respec- 
tively. The Matrix B is defined as 



1 

-1 a 



N-l 



(6) 



Here, N is the number of unit cells, c.f. Fig. 1. The 
variables a, /?2 and a are introduced as 



Ae 



e + Ae + Tir 2 



Ae 2 



ft 



e + Ae 
e 2 - Ae 2 - 



nT 2 
^2 



(7) 
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where e represents the energy of the incoming particles 
and is defined in units of ts identically to e{n) in (2). The 
wide-band-limit (energy independent lead properties) is 
used here for the <S-matrix. The matrix-elements of B 
essentially control whether the lattice is a conductor or an 
insulator, because it is the quantity that mainly controls 
the denominator of S12 (and therefore the transmission 
probability T = | S12 1 2 ) and contains the TV-dependency. 
The parameter a determines whether the elements of B 
grow exponentially with N or oscillate. If the definition 
of a in (7) is solved for e, the spectrum of the infinite 
lattice (2b) is recast with a replacing 2 cos 



ei,2 



a 



± 



Ae 



• Tir 2 a . 



(8) 



One can deduce that if and only if the value of e is part 
of the energy spectrum, there is a real valued n satisfying 
a = 2 cos ft and hence \a\ < 2. An elementary analysis 
of matrix B yields that its elements grow exponentially 
with N if I a | > 2. With \a\ < 2, the dependency is 
periodic and at \a\ = 2 it is linear. Respective explicit 
expressions for B are given in Eqs. (A15) and (A16) in 
the Appendix. We conclude that transmission T = | S12 1 2 
through the finite lattice is suppressed exponentially with 
N if the particle energy is not part of the spectrum of the 
infinite lattice. For energies that belong to the spectrum 
the transmission oscillates and the number of resonances 
is related to N. The transmission is visualized together 
with the parameter a as a function of the particle energy 
e in Fig. 3. For one example of energy bands [Fig. 3(a)] 
the transmission is presented for N = 3 [Fig. 3(b)] and 
N = 9 [Fig. 3(c)]. The parameter a, which mainly con- 
trols the transmission, is shown in Fig. 3(d). One can 
see that the interval a G [—2,2] corresponds exactly to 



the spectrum. Transmission outside of this interval is 
suppressed as a function of N [compare Fig. 3(b) and 
3(c)]. Figure 3(h) illustrates the effect on the transmis- 
sion when a gap opens in the energy spectrum, namely 
the suppression of the transmission inside the gap. 



C. Adiabatic quantum pumping properties 

Having the scattering matrix of the finite triangle lat- 
tice, the scattering approach to the AQP 2 can be ap- 
plied. By varying slowly two of the three parameters n, 
T2 and Ae, a net current can be produced in absence of 
an external bias. In the adiabatic regime the current is 
proportional to the frequency uj of the periodic variation. 
The adiabatic approach is exact in the limit uo — > 0. We 
calculate the charge transfer Q, which is independent of 
(j. It is the integral of the current over one parameter 
cycle. Brouwer's formula for the charge transfer in the 
zero temperature limit reads 2 



Q 



dXidXj 3 < Tr 



dS <9S f 

~dX~ 2 dX[ 



. (9) 



where j is the lead-index and X\ and X2 are the pumping 
parameters. The trace is applied over the channels and 
can be omitted in the case of one-channel-leads. The 
scattering matrices have to be evaluated at the Fermi- 
energy, which will be expressed with the dimensionless 
scaled parameter ep in analogy to the scaling of e{n) 
in (2). The choice of the lead j determines the sign of 
the charge transfer. It is convenient here to choose j = R 
which fixes the result to the transfer from the left to the 
right side in Fig. 1. 

There are three possible choices of pumping parame- 
ters: {ti,T2}, {ri,Ae} and {r2,Ae}. Exchanging the pa- 
rameters r\ and T2 is equivalent to the exchange of the 
orientation along the lattice axes. Hence, it is equivalent 
to a sign change of the charge transfer. It also follows 
that the charge transfer vanishes if t\ = T2 holds for the 
whole cycle. After fixing the set of pumping parameters, 
we study how the choice of the parameter cycle influences 
the charge transfer. In the following we will focus on the 
choices: Path A =>• {ri^} and Path B =^> {r^Ae}. 

For both cases the charge transfer is plotted as a func- 
tion of the extreme of the pumping parameter path and is 
always expressed in units of the electron charge (e = 1). 

Figure 4 shows examples where parameters Path A is 
varied around band touching points [c.f. Fig. 2(a) and 
3(a)]. Here the lower panels show the transmission prob- 
ability together with the pumping parameter cycles. In 
the upper panels the charge transfer is plotted versus the 
maximum value of T\ in the cycle. In the following we 
summarize the main features of the of the transferred 
charge Q along paths A and B: 

1. Crossing of a spectral gap — If the pumping path 
ends in a gap of the finite system energy spectrum, 
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FIG. 4: (Color on-line) Varying parameters n and Ae — path B — around two band touching points (marked with black 
crosses) for different Fermi-energies cf and N. The lower panels show the parameter cycle (blue) on the density plot of 
transmission T (black: T = 1, white: T = 0). In the upper panels the charge transfer is plotted versus T ( max ). Vertical dashed 
lines (cyan) indicate the boundary between insulating and transmitting regions. For all panels: T2 = 8 



the charge transfer Q is quantized to an integer 
value as shown in Figs. 4(a)-4(c) and 5(b). 2,13 ' 19 
If the pumping path end in one "zero-gap" of the 
spectrum, thus a Dirac-like points of the energy 
spectrum [crosses in the lower panels of Fig. 4], the 
charge transfer Q is also quantized but to an half- 
integer values so as observed also in graphene. 20,21 

2. Dependance on the number on unit cells N — The 
charge transfer Q can either be an increasing or 
a decreasing function of the number of unit cells 
N. The exact dependence depends on the system 
parameters. 22 Figure 5(c) shows for path B travers- 
ing configurations with nonzero transmission. The 
charge transfer for this cycle is calculated for differ- 
ent numbers of iV, which effectively represents the 
length of the lattice, and is visualized in Fig. 5(d). 
On average, there is a linear increase of the charge 
transfer with respect to N. 

3. Influence of the coupling to the leads c — Both in 
the case of very weak or very strong coupling c the 
pumped charge is quantized to a constant value. 
This is because in both cases the parameter paths 
are ending in the energy gap of the of the finite size 
system. 



4. Crossing of a flat band — In the case of path B it is 
possible to cross a flat band of the infinite lattice, 23 
[c.f. Eqs. 3]. In Fig. 5(a) and 5(b) the parameters 
are varied around the fiat band. The bands cross at 
n = r 2 = if |Ae| < 1 [see Fig. 3(e)]. This corre- 
sponds to the B-site (flat band) completely decou- 
pled from the rest of the chain (cosine band). How- 
ever, if the parameter cycle surrounds this point 
[see Fig. 5(a)], the structure of the whole system — 
with nonzero coupling r± and T2, e.g. of Fig. 3(g) 



— is then responsible for the charge transfer. Fig- 
ure 5(b) shows that the charge transfer becomes 
Q = — 2e for large pumping cycles that purely tra- 
verse parameter configurations where transmission 
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FIG. 5: Demonstration of different charge transfer results 
with pumping parameters n and T2 — path A. Panels (a) 
and (b): Symmetric pumping cycle with r^ max ^ = — r^ min \ 
Ae = 0.5, €f = and N = 11 around crossing point. Pan- 
els (c) and (d): Charge transfer of pumping in transmitting 
region increases linearly with lattice length N. Ae = 1.8, 
cf = 1. Brightness of panels (a) and (c) corresponds to trans- 
mission T (Black: T = 1, White: T = 0). 
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is inhibited. The potential benefits of crossing a 
flat band 23 is always cancelled by the fact that a 
parameter path of type A is always crossing this 
flat band two times. 



III. CONCLUSION 

We have presented the general spectral properties of 
the triangular lattice. We have shown that by chang- 
ing the system parameters, the spectrum show one- 
dimensional Dirac cones, flat band and many other fea- 
tures. The transport properties of the finite length lattice 
are also largely affected by the system parameters. 

Thus, within the scattering matrix approach to the adi- 
abatic pumping a large set of results can be produced. It 
turns out that, if the parameter configurations of the cy- 
cle remains in an insulating area (no transmission), the 
charge transfer is quantized — either to an integer or 
an half-integer value. At zero transmission the Fermi- 
energy is positioned between the bands (without touch- 
ing them). In this case, the scattering matrix only con- 
tains nonzero reflection amplitudes that vary along the 
parameter cycle. Although all particles are reflected in 
the static case, the variation in time allows for scattering 
processes that lead to a net charge transfer Q through the 
system that is independent of the length of the lattice. 

If the parameter path crosses areas with nonzero trans- 
mission, the charge transfer is generally a non-integer 
value. If transmission is allowed, the charge transfer de- 
pends sensitively on all used parameters. If particles are 
transmitted through the lattice, the amount energy - 
gained or lost — depends on TV, the number of unit-cells 
of the lattice [c.f. Fig. 5(d)]. This is related to the pa- 
rameter change being applied to all parts of the lattice. 
A perturbation affecting only one specific unit cell of the 
lattice does not lead to an TV-dependency. 

Note that this linear dependency on TV does not mean 
the generation of an arbitrary large current with increas- 
ing lattice size. Because of the approximations used in 
this model the pumping frequency uj needs to be suffi- 
ciently small. This requires uj to decrease with increasing 
length of the quasi-one-dimensional lattice. 
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Appendix A: Calculation of the S-matrix 

The Hamiltonian in position space has the form 

/H L T, 0\ 
H = Tl H S tH , (Al) 
V T 2 H L / 

where Hl is the Hamiltonian of the (semi-infinite) linear 
chain and H$ is the one of the lattice sample. Ti and T2 
are the matrices that connect lead sites with sample sites 
and here have only one non-zero element. The retarded 
Green's function, defined with 

[(£ + i77)l-H]G R (£) = l , (A2) 

can also be partitioned into sub-matrices 

(Gn Gis Gi 2 \ 
Gsi G s G S2 • (A3) 
G21 G2S G22 / 

The Fisher-Lee-relation relates the Green's function ele- 
ments to the elements of the scattering matrix. 16 18,24 

The leads are modeled as linear chains with dispersion 
relation E K = 2£ 3 cos k + e^. It is useful to consider the 
wide-band-limit, where k is fixed and the Green's func- 
tion does not depend on E K or on on-site energy e^. If 
£3 ^> E K — €l, the choice k = | (cos7r/2 = 0) is valid 
for all particle energies. The Fisher-Lee-relation for this 
case can be expressed with 

S = -i + 2it 3 ( £s)N lett>Nlett (G s )jv left) iv right \ _ (A4) 

ight bright / 

The elements of the Green's function are restricted to the 
translational invariant region. Here, the edge-sites of the 
triangle lattice can be considered as the beginning of the 
translational invariant leads. Four elements of Green's 
function Gs are needed to calculate the scattering ma- 
trix. To determine these quantities it is useful to extract 
from (A2) three sub-matrix equations, 

[(E + i77)l-H L ]Gis-TiG s =0 

[{E + irj)t - H L ]G 2S - T 2 G S = (A5) 

[(E + ir])t - H S ]G S - TlGis - TjG 2S = 1 

They can be used to derive 

G s = [(E + irj)t - H S - Tlg L Ti - T^Ta]" 1 

(A6) 

= [(E + i V )l-H eS ]-\ 

where gL = [{E + iry)! — Hl] -1 is the Green's function 
for a single isolated lead. Its edge element in position 
space reads 15 

(kl)h = r G ~ iK(E) = -r> (A7) 
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where the second identity holds within the wide band 
limit. The effective Hamiltonian in position space reads 



(Heff)*j = (H s )»j - (&iN left 8j Nleft + S iN . ht S. 



ght j bright 



)i* 3 



To demonstrate its symmetry and in order to provide an 
appropriate calculation of G$ the effective Hamiltonian 
is divided into sub-matrices as 



H e ff = 



/Hl 


T 





... o 


^ 


T t 


H c 


T 









T t 


H c 


'•. 











. T 













T f H c 


T 









T f 


Hr/ 



(A8) 



with 



H L 

Hr 



ti 6b 

e A - it 3 




Hc = 
T = 



h o 
t 2 o 



Now we can express the relation 

[(£ + i77)l-H eff ] G s = l 

as 



(A9) 



[(E- 


Mj7)1 


-H L ] 




- T G 2j 


= 


PH 


-ir?)l 


-H c ] 










-ir?)l 


-Hr] 


G^ 


- T^Gjv-y 





where z e [2,7V - 1] and j e [l,iV]. Equation (A9) 
contains twelve equations for the elements of the sub- 
matrices which can be solved independently for every in- 
dex j. Only the sub-matrices Gn, Gin, Gni and Gat at 
are relevant for the scattering matrix. And only the (1,1)- 
element (A-site) is needed from the sub-matrices. By in- 
troducing gi = ts(Gij)u, the set of relevant algebraic 



equations can be simplified to 

Pi9i ~92 = Sjia (Alia) 

&9i ~ 9i+i ~ 9i-i = (Allb) 

fcgN - 9n-i = S jN a, (Allc) 

where a, a, f3\ and fa are defined as in Eq. (7). Further, 
we express Eq. (Allb) as 



9i+i) = ( 1 
g i+2 J \-l a) \g i+ il ' 



(A12) 



which allows to connect the last two elements with the 
first ones as 



9N-1 ] = B (9i 
9n 



with B = ( °, 1 
92) ' V-l ol 



N-l 



From Eqs. (All) and (A13) it follows that 
5ji (B12 — AB22) — SjN 



9i=cr 
9n = o- 



Bii + (/3i+/3 2 )B 12 -/3i/3 2 B 22 

SjN (B12 - feB 22 ) — 5ji 
Bn + (^1+^12-^1^22' 



(A13) 

(A14a) 
(A14b) 



which yields the scattering matrix ((5)). 

The matrix B can be evaluated analytically to read 



1 f-an((N-2)u) sin((AT - 1)^ 



with cos v = a/2 for the case \a\ < 2. On the other hand, 
for \a\ > 2 we obtain 

_ _JL_ (- sinh((7V - 2)1/) sinh((7V - 1»\ Am) 
~ sinh v \~ sinh((A^ - \)v) sinh(AT^) y 

with coshz/ = |a/2|. 
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